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ABSTRACT 



We report a mass and rotational broadening (vsini) for the pulsating white 
dwarf component of the WZ Sge type Dwarf Nova GW Lib based on high- 
resolution VLT spectroscopy that resolves the Mg II 448lA absorption feature. Its 
gravitational redshift combined with white dwarf mass-radius models, provides 
us with a direct measurement of the white dwarf mass of Mi = 0.84 ± O.O2M . 
The line is clearly resolved and if associated with rotational broadening gives 
vsini = 87.0 ± 3.4 km s -1 , equivalent to a spin period of 97 ± 12 s. 

Subject headings: binaries: spectroscopic — novae, cataclysmic variables — stars: 
individual (GW Lib) 
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Introduction 



The population of cataclysmic variables (CVs) provides an important sample of binary 
systems possessing homogeneous configurations of white dwarfs (WDs) accreting from near 
main-sequence donor stars. A CV initially evolves towards shorter orbital periods, but near 
a period of ~ 75min bounces back towards longer periods with its donor star turning into a 
degene rate brown dwarf. Recent CV searches, s uch as the significant harvest of CVs from 



SDSS ( Szkody et al. 



2009 



Gansicke et al. 



20091 ) are finally unearthing large numbers of 



systems at short orbital periods as expected from binary evolution considerations. 

Reliable estimates for binary parameters are needed to place individual systems on 
their evolutionary tracks. The faintness of the low mass donor stars in comparison to the 
WD and its a ccretion flow makes s uch parameter estimates difficult for short period systems 



(e.g WZ Sge; 



Steeghs et al. 



2001 



Steeghs et al 



eclipse constraints can be exploited (ILittlefair et al. 



2007) , ex c ept in favourable cases where 



200J). 



Patterson 



(120011 ) has identified 



an indirect method of inferring binary mass ratios using the tidally driven superhump 
modulations in the lightcurves of CVs. Although promising, this furthermore highlights 
the need for suitable calibrator systems for which binary parameters can be accurately 
determined independently. 



GW Lib is a short period dwarf nova (P or b=76.78 mins; iThorstensen et al.l 120021 ) 



displaying the typical characteristics of a CV near its p eriod minimum. C rucially, it 



2004 ). by analogy 



was the first CV discovered to contain a pulsating WD (Ivan Zyl et al. 
with field WD pulsators that lie in the instability strip for non-radial pulsations. So far, 
various methods to determine its system parameters have been expl oited, including mo del 



fits t o the WD abs orption profiles, pro minent in both the optical (IThorstensen et al. 



20021 ) and the UV (ISzkodv et al. 



20021 ) providing estimates for the WD T e g and log p. 



Asteroseismological models, combined with a UV-flux limit, suggests a WD mass of 
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Mi = L02Mq ( 



2007 flTemp 



Katoetal 



2005|). 



cton 



bwnslev et al.l 120041 ) . After its second recorded super-outburst in April 



20071 ). the superhump period suggest a mass ratio of q = M?jM-\ = 0. 06 



20081 ) when combined with Patterson's empirical relation (IPatterson et al. 



In this letter we explor e the possibilities of u s ing th e gravitational redshift of the Mg II 

( 201(3 ) in combination with mass-radius 



absorption line reported in 



van Spaandonk et al 



relations to give an independent measurement of the mass and spin of the WD in GW Lib. 



2. Observations and reduction 

GW Lib was observed within the 69.D-0591 program at the 8.4-m Very Large 
Telescope (vlt), located at the Paranal Observatory in Chile as part of the European 
Southern Observatory, equip ped with only the blue arm of the Ultraviolet and Visual 



Echelle Spectrograph (uves: 



Dekker et al. 



20001 ) . We retrieved a total of 45 science frames 



obtained during 2002 May 16-17 covering 2.4 and 2.8 binary orbits on the respective 
nights. The reduction of the raw frames was conducted using the most recent standard 
recipe pipeline release of the UVES Common Pipeline Library (cpl) recipes. The resultant 
optimally-extracted spectra covered a wavelength range of A4020-5240A at a dispersion 
of 0.031A pixel -1 and a spectral resolution of 0.10A (5.77 km s -1 ) as measured from the 
skylines. The spectra were wavelength calibrated with one Thorium Argon arc per night. 
This calibration was tested against the two sky lines visible at A5197.92A and A5200.28A. 
The corresponding science frames were corrected for any remaining shifts. The residuals 
were scattered around zero with a maximum amplitude of 0.01 km s _1 . Next, heliocentric 
velocity corrections were applied to the individual frames to deliver spectra in a common 
heliocentric rest frame. 
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No standard was observed in the correct settings on the nights and as no master 
response curve exists for the non-standard setup used, the frames were not flux calibrated. 
The exposure time was 500 seconds giving a signal to noise of ~6.5 per spectrum. Details 
of the observations can be found in Table [1] while Figure [1] shows the average spectrum of 
GW Lib on the 16th of May. Prominent features are the Balmer disc emission lines on top 
of broad absorption troughs from the WD, visible due to the low mass accretion rate in the 
system. He I is also seen in emission as is He II at 4685. 75A. 

Our pre-outburst, inter mediate resolution spectra o f GW Lib showed Mg II at 4481. 21 A 



in absorption (see figure 2, 



van Spaandonk et al 



2010). The archival vlt/uves data 



confirm the presence of this line and thanks to the superior spectral resolution combined 
with the low inclination of the system shows it to be unblended from the nearby He I 
emission at 447lA. We measure an EW of 0.25 ± O.OlA (similar to 0.24 ± 0.03A in 2004) 
and a Full Width at Half Maximum (fwhm) of 1.35 ± 0.04A. 



3. Gravitational Redshift for Mg II 



In low accretion rate dwarf novae, the luminosity of the accretion disc is low enough to 
show the broad absorption features of the WD and can even show the narrow absorption 
features of metal lines due to freshly accreted gas. These lines open a window to probe the 
WD atmosphere directly and give independent measurements of stellar parameters. For 
lines formed near the primary, a gr avitational red-shift is expected, introduced in the deep 



gravi tational potential of the WD ( lEddington 



1924 



Greenstein fc Trimble 



1967 



Sion et al. 



19941 ) . A measurement of the gravitational redshift in the rest frame of the binary could 
provide the WD mas s directly when combined w ith mass-radius models (e.g. Eggleton's 



relation as quoted in 



Verbunt &: Rappaport 



in CVs in the cases of U Gem (jLong fc Gilliland 



1988 



This method h as previously been used 



19991 ). VW Hyi f lSmith et al 



200J) and 
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WZ Sge flSteeghs et all 120071 ) 



In the case of GW Lib, the directly measured redshift of the magnesium line (vMgii) 
needs to be rectified for several contributions in order to get the true gravitational redshift 
induced by the WD only (% av (WD)). These corrections consist of the systemic velocity 
(7) of the binary system and the effects of the gravitational potential from the donor star 
(fgrav (donor)). Hence the gravitational redshift due to the WD is given by: 

fgrav(WD) = v MgU - w grav (donor) - 7 



We will discuss the various contributions and measurements independently. Firstly, all 
45 spectra of the two nights were individually continuum normalised and then binned into 
20 equally spaced orbital phase bins to increase S/N. 



3.1. t>MgII 

To measure the redshift of the magnesium line to the best precision and minimise any 
orbital effect we compared several methods. 

Firstly, we made a weighted orbital average and fitted the Mg II absorption line with a 
weighted triple Gaussian as the Mg II is a triplet line. The different components have rest 
wavelengths of A4481.126A with a transition probability of log(of) = 0.7367, A4481.150A 
with hg(gf) = -0.5643 and A4481.325A with hg(gf) = 0.5818lJ. The variables for the 
fit are the common offset (for all 3 lines) and the common fwhm to provide a good fit to 
the blue wing of the absorption feature. The peak height is also a common variable but 



x Data from The Atomic Line List Version 2.04: 

http : / /www . pa . uky . edu/~peter/ atomic/ 
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scaled according to the various transition probabilities. This gives a best fit with a mean 
offset of 35.2 ±1.1 km s _1 , a fwhm of 1.32 ± 0.08A and a common absorption line depth 
of 2 ±0.1%. 

Secondly, we fixed the fwhm and peak of the Gaussians to these values and checked 
the individual spectra for orbital variablity. Unfortunately, the S/N is too low to give good 



individual fits (individual la errors on the offset are ~ 8 km s 1 ) nor can we phase- 



ock this 



van Spaandonk et al 



fcOloh due to 



motion to the ephemeris of the system as determined in 
the uncertainty in the period. The resultant radial velocity curve suggested motion, and a 
formal sine fit delivered a semi-amplitude of K = 13 ± 2 km s _1 . To rectify for any orbital 
motion, we removed any measured shift compared to the mean from individual spectra. 

An uniform orbital average was constructed to minimise any effects caused by varying 
S/N over the orbital period and finally we refit this last average with our triple Gaussian 
fit with as variables the common offset, the common fwhm and the weighted peak as 
described before (see inset Figure [Q . 

This gives a final measurement of t>Mgii = 35.8 ± 1.5 km s _1 . The line has a fwhm 
= 1.30 ± 0.05 A. We have tested the accuracy of our measurement by following different 
recipes for combining and averaging the spectra but note that the uncertainty on the 
gravitational redshift is dominated by the S/N and resolution of the data and not the 
specific recipe used. 



3.2. f grav (donor) 

The first correction is due to the influence of the gravitational potential of the donor on 
the magnesium line. However, the expected low mass of the donor star (~ 0.05 M Q ) gives 
only a small correction of 0.06 ± 0.02 km s _1 near the WD surface, effectively negligible 
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given the measurement uncertainty. 



3.3. The systemic velocity 7 



From simultaneous double Gaussian fits to the double peaked H/5 and H7 disc lines, 
combined with radial velocity curve fits provide a systemic velocity of —12.3 ± 1.2 km s _1 



and a semi-amplitude of Kd^r. = 36.4 ± 1.8 



•cm s . These values are consistent with 



previously derived values (IThorstensen et al 



2002 



van Spaandonk et al 



2010|) 



As we measure 7 from disc lines, we need to take into account that this emission is 
red-shifted both by the WD and the donor star. As t> grav is reciprocal to the distance, the 
amount we have t o account for i s minimal at the edge of the disc, -Router disc — 2.2 x 10 8 m 



(Equation 2.61 of 



Warner 



19951 ) giving a minimal correction to v grav of 1.7 km s 1 . 



A more realistic value comes from the projected Keplerian velocity at the edge of 



the disc K^ sc ~ 200 km s 1 (from the location of t 



and the estimate for i, see 



van Spaandonk et al. 



re disc ring in Doppler maps 



2010) which gives a Kepler speed of 



■^disc = K&\scl srn ^ ~ 1070 km s _1 and corresponds to a t> grav = v 2 /c ~ 3.8 km s _1 . Including 
the transverse Doppler redshift (at half the strenght of the gravitational redshift) this 
amounts to a total gravitational redshift at the location of the disc lines of 5.7 ± 1.6 km s _1 . 
The gravitational potential of the donor star has an effect of only 0.08 ± 0.02 km s _1 at a 
distance of ~ 3 x 10 s m. Hence the total t> grav (disc) = 5.8 ± 1.6 km s^ 1 . 

The 7 from the disc lines combined with the above correction should be consistent 
with the 7 suggested by the radial velocity curve from the donor star (corrected for 
the effects of the gravitational potential at its surface). F rom the Ca II emissio n 
line in the i-band we previously found 7 = —13.1 ±1.2 ( Ivan Spaandonk et al. 



20101 ). the correction at the surface of the donor star is 0.34 ± 0.15 km s 1 from 
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the donor star and 1.4 ± 0.2 km s 1 from the WD (again including the transverse 
Doppler shift). Thus we have 7 disc = (-12.3 ± 1.2) - (5.8 ± 1.6) = -18.1 ± 2.0 and 
7donor = (—13.1 ± 1.2) — (1.7 ± 0.3) = —14.8 ±1.2. Both estimates for 7 are indeed consistent 
and with similar precision. Hence we have used 7di SC as it is derived from the same data set. 

3.4. The implied WD mass 

Combining all values we find for the final gravitational redshift: 

%av(WD) = (35.8 ± 1.5) - (0.06 ±0.02) - (-18.1 ±2.0) kms" 1 
= 53.8 ±2.5 kms" 1 



When combined with theoretical and empirical models for the mass-radius relationship for 



WDs this gives a direct measurement of the WD ma ss. The models we used are 
zero-temperature mass-radius relation as quoted by 



Verbunt fc Rappaport 



Iggletons 



(1988) and 



Fontaine et al 



(2001) 



several appropriate non-zero temperature models for GW Lib from 
We plot these models in Figure [2] together with the M 1 (R) line demanded by our measured 
f grav . Accommodating the intersections between non-zero temperature models we find that 
Mi = 0.84 ±0.02 M Q . 



3.5. System parameters and WD spin 



Kato et al. 



( 120081 ) reported the detection of superhump modulations in GW Lib. 



Combining this with the period from lThorstensen et al. 



excess - mass ratio relation given by 



( 120021 1 and the improved superhump 



Kniggd ( 120061 ). the implied mass ratio of the system 



is q = 0.06 ± 0.008. We previou s ly det ermined the projected radial velocity of the donor 
star, K 2 in Ivan Spaandonk et al.l J2OI0I ) . Thus with an independent determination of M\ 
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from this paper, we can solve the system parameters for G W Lib under these constraints 
and list these in Table [2J 

If one assumes that the width of the Mg II absorption line is dominated by rotational 
broadening, its fwhm can constrain the spin period of the WD. The fwhm of the absorption 
line fit is measured to be 1.30 ± 0.05A which corresponds to a vsini of 87.0 ± 3.4 km s -1 
at this wavelength. Assuming i is close to the value giving in Table |2l this translates 
into a rotation speed of 448 ± 24 km s _1 at the surface of the WD. For a WD of mass 
0.84 ± 0.02 M & and a radius for 6.95 ± 0.15 x 10 8 cm this results in a spin period of the 
WD of 97 ± 12 seconds. 



4. Discussion 

We have measured the gravitational redshift of the Mg II absorption line in high- 
resolution echelle spectra of GW Lib during quiescence. Assuming an origin in the 
photosphere of the accreting WD, we combined this redshift with non-zero temperature 
mass-radius relations for WDs to derive a WD mass of 0.84 ± O.O2M . Combining this 
independent measurement with other constrains confirms that GW Lib is a low mass ratio 
system observed at very low inclination (Table [2]). Because the Mg II line was well- resolved 
in our data, we could also estimate the spin period of the WD to be 97 ± 12 seconds if we 
assume the width of the line to be dominated by rotational broadening. 



Our mass value is significantly below the 1.02 M & lower limit derived by 



Townsley et al. 



( 12004] ) . However, their choice of prefe rred solution was lar gely driven by the WD size 



as implied by the measured UV flux. 



Szkody et al. 



(120021 ) fitted WD models to UV 



spectra of GW Lib. Two different single temperature solutions were explored; d = 171 pc 
(Mi = O.6M /v grav ~ 29 km s" 1 ) and d = 148 pc {M 1 = O.8M /w grav ~ 49 km s" 1 ) 
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(Thorstensen, private communication - improvement of the distance in 



respectively. These solutions can be corrected to the latest parallax dis tance of lOOl }, pc 



Thorstensen et al. 



20021 ) . At this distance, the observed UV flux implies a radius of 5 x 10 8 cm (for a WD 



temperature of 14 700 K). However, such a small radius intercepts th e mass-radius re l ations 



at 1.07 Mq, Figure EI (horizontal solid line) and was indeed the reason 



Townsley et al. 



(120041 ) 



ruled out their lower mass solutions. Such a massive WD would result i n t>g. rav ~ 



95 km s 



Szkody et al. 



fl2002h 



much larger than our measured redshift would suggest. We note that 
fail to find a good single temperature fit and claim the best fit arises fitting a dual 
temperature model to the spectrum, allowing for a cooler zone that would place GW 
Lib within the ZZ Ceti instability strip for pulsating single WDs. However, the UV-flux 
versus implied WD radius then becomes a less clear-cut argument given the freedom of two 
temperatures and the surface coverage split between the two. We note that the reverse 
has been seen in U Gem wher e the UV-flux places the WD at a lower mass than the 



gravitational redshift suggests (jLong fc Gilliland 



1999 



Long et al 



20061 ) 



As 



Townsley et al. 



(120041 ) warn, no WD rotation was included in their asteroseismology 



models and therefore their derived parameters (Mi, M acc , M) may need to be revisited given 



our estimate of the spin period of th e WD in GW Lib. A P< 



give rise to rotationally split modes. 



van Zyl et al. 



spin 



97 ± 12 seconds should 



(120041 ) conducted a thorough campaign 



with a baseline of over 4 years to find these modes in the power spectrum of GW Lib. 
The extensive search revealed a possible doublet around the 230s mode with a frequency 
difference between the components of 0.79 /xHz. If this doublet originates from the rotation 
splitting of an I = 1 mode, the WD would have a spin period of ~ 7.3 days which is much 
longer than expected for the WD in a CV. Approached from the other side, if a similar 
mode is split as a result of a rotation period of ~ 100 s, the frequency difference would 
be ~ 0.5 mHz and should be visible in the power spectrum. Note that for the proposed 
spin period, the first order approximation for frequency splitting is no longer valid and the 
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splitting would become asymmetric. 

From both GW Lib's and U Gem's comparison between WD mass found by 
gravitational redshift and the constraints from the UV-flux we can conclude that the 
measurement of the gravitational redshift through WD absorption lines gives in principle a 
good measurement of the WD mass in CVs, but is not always fully consistent with other 
constraints. Nonetheless, it suggests photospheric absorption lines in the atmospheres of 
accreting WDs can be a viable tool for mass measurements, in particular for low mass 
accretion rate CVs where the light from the WD dominates and especially for low % systems. 

For GW Lib itself, we are still dependent on a number of assumptions in order to 
solve its system parameters. To properly align all methods, we need to revisit the mass 
from asteroseismology taking into account the effects of rotation and see if a WD model 
with a mass of O.84M at the latest distance can be fitted to the UV data with reasonable 
temperatures. To be able to calculate a direct mass ratio, a determination of K 1 from the 
Mg II should be possible through high S/N phase- resolved spectroscopy. If furthermore 
matched with an solid ephemeris using the emission from the donor star (giving K 2 and 
7), a fully consistent and accurate WD mass for this prototypical accreting WD pulsator is 
entirely feasible. 

DS acknowledges an STFC Advanced Fellowship. TRM was supported by an 
STFC Rolling Grant. Based on observations made with ESO Telescopes at the Paranal 
Observatory under programme ID 69.D-0591 
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Table 1: Details of observations 



Date 


n spectra 


Exp time 


Start/Finish 


Orbital Coverage 






(s) 


UTC 


(^obs/ -forb) 


16 May 2002 


20 


500 


05:28:26 - 08:34:02 


2.4 


17 May 2002 


25 


500 


04:49:03 - 08:24:08 


2.8 



- 15 - 



- 0.9 

- 0.8 



He I 
I 

I 



4100 4200 4300 4400 4500 4600 4700 4800 4900 5000 5100 5200 

Wavelength A 



Fig. 1. — Average spectrum of GW Lib on the 16th of May 2002. The inset shows the Mg II 
absorption line averaged over 1 orbital period with the final Gaussian fit. The total fit (solid 
line) is the sum of three separate fits to the triplet lines, all according to their own transition 
probability. A4481.126 dotted, A4481.150 dotted- dashed, A4481.325 dashed. 



Table 2: System parameters for GW Lib 



P (min) 


76.78± 0.03 


Thorstensen et al. 2002 


Mi (M ) 


0.84 ±0.02 


This paper 


Q 


0.060 ±0.008 


Kato et al. 2008: Knieee 2006 


K 2 (kms- 1 ) 


100.8 ±7.1 


van Spaandonk et al. 2010 


K x (kms- 1 ) 


6.25 ±0.4 


From q and K 2 


M 2 (M ) 


0.050 ±0.007 


From q and Mi 


Vi (kms -1 ) 


30.8 ±0.5 


Kepler 3rd Law 


t>2 (kms" 1 ) 


513.4 ±4.8 


Kepler 3rd Law 




11.2 ±0.4 


From K 2 and v 2 


-fspin(s) 


97 ± 12 


From v sin % 
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0.4 - 



°t?.75 0.80 0.85 0.90 0.95 1.00 1.05 1.10 
M(M ) 

Fig. 2. — Mass radius relation ship for the measured t> grav = 53.8 km s" 1 (dashed line 
from lower left to upper right) with la error lines (dotted lines). The mass-radius models 
are Eggleton's zero-T relation (solid line from upper left to lower right), a T = 15 500 K, 
qH = 10" 4 and gHe = 10~ 2 C-core model (dashed line). Models with the same core but 
diferent temperatures at T = 18 000 K and T = 13 000 K are similar and hence not plotted. 
For a WD with a CO-core, T = 15 500 K, qR = 10~ 10 and gHe = 10" 2 , the model is almost 
similar to the zero-T relation. The horizontal lines denote the UV flux estimate from the 
Szkodys paper corrected for the latest distance measurement for GW Lib (solid line) and la 
errors (dotted) 



